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I. INTRODUCTION 



One of the fundamental mysteries of neutrino-physics is the nature of massive neutrinos 
They can have Dirac or Majorana masses [^J, with their own signatures in each case. 
In particular, it is assumed that the best way to establish their character is to search for 



neutrinoless double f3 decay ((PP)q u ) 



oest way 

0, Q, fl. 



Indeed, its mere observation would imply 



that neutrinos are Majorana particles. However, if neutrino masses are a fraction of eV and 



no New Physics (NP) gives large enough contributions, it will be quite difficult to observe 
{P/3)ov flflflS- 



Although oscillation experiments are not sensitive to the absolute scale of neutrino masses, 
they are able to probe very tiny neutrino mass differences Q, ll|, 12]. Hence, the obvious 
question is if neutrino oscillations can distinguish between Dirac and Majorana neutrino 
masses . This is not possible in neutrino oscillations in vacuum [15l. Ilq|. The reason 

is that the additional phases in the neutrino mixing matrix distinguishing Majorana from 
Dirac neutrinos disappear from the oscillation probability expressions. The oscillation in 
matter can differ from that in vacuum, leading to both the change of the effective mixing 
angles and of the effective masses 17|, ]l8|, ]l9|- However, if neutrinos are relativistic and only 



interact through left-handed (LH) currents, the oscillation in matter does not discriminate 

Detween both types of fermions either. This is not the case for non-relativistic neutrinos 
1 

201 ] but unfortunately these neutrinos are not available in any experiment, at least up to 
now. So, within the New Standard Model (vSM) 21j and for neutrinos in the mass range 
of tens of meV it will be very difficult to find experimental evidence which can distinguish 
Dirac from Majorana neutrinos. 

Then, one may wonder if there is some type of NP which can modify neutrino propagation 
in a dense medium 22|, ]23|, ]24], ]25|, ]26|] to allow for the determination of the neutrino mass 



character. As we will review, although NP can give rise to a low energy Hamiltonian with 
any tensor structure for neutrino bilinears uT a u, with T a = I, 75, 7 M , 7^75, and a^, the 
scalar and pseudoscalar currents vanish or are suppressed for relativistic neutrinos. They 
are O(^), where m is a light neutrino mass and E the neutrino energy. Vector and axial- 
vector helicity flip and tensor helicity non-flip transitions have similar suppression factors. 
Only the vector and axial-vector helicity non-flip transitions and the tensor helicity flip 
ones are unsuppressed, and may then help to discriminate between relativistic Dirac and 
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Majorana neutrinos propagating in a dense medium. Here we will concentrate on the spin 
non-flip vector and axial-vector transitions. Tensor interactions have been discussed in this 
context in [^tJ. In any case, as the production and detection of neutrinos (antineutrinos) 
with positive (negative) helicity are strongly suppressed, the sensitivity to spin-flip tensor 
transitions is reduced. Besides, they require a polarised background to manifest, and if 
generated at higher orders in perturbation theory, they are also extra suppressed by loop 

, n 

factors 1 [28|]. 

In this paper we will consider an effective Lagrangian with only arbitrary vector and 
axial-vector couplings modifying the LH structure of the vSM. This is enough to generate 
(new) vector and axial-vector four-fermion interactions which are those to be probed by 
relativistic neutrinos. The new couplings affect differently the propagation of Dirac and 
Majorana neutrinos in matter. The largest effect results from the addition of new right- 
handed (RH) neutral currents, being the corresponding transition probability dependence 
linear 2 . We discuss its size for different channels as a function of the neutrino energy and 
the baseline, showing that it can eventually allow for discriminating between Dirac and 
Majorana neutrinos. 

We first revise the necessary known results on neutrino propagation in a dense medium. 
In next Section we describe the neutrino interaction with matter introducing a general, 
relativistic invariant Hamiltonian, and show that for relativistic neutrinos only vector, axial- 
vector and tensor terms have contributions to neutrino oscillation amplitudes which are not 
suppressed by ^ factors. In Section (Till we discuss the Lagrangian with general LH and RH 
currents. As deviations from the vSM are strongly constrained by existing experimental 
data, they must be small. Then, using the corresponding four-fermion effective Hamiltonian 
we calculate the effective interaction with matter for Dirac and Majorana neutrinos. In 
Section IIVI the probabilities for neutrino oscillations in matter are derived for a neutral 



background with uniform spin and momentum distributions 



291 ] . We present our numerical 



results for neutrino propagation inside Earth for different channels, baselines and neutrino 



1 They also result from Fierz rearrangement of (pseudo) scalar interactions and can be of the same order as 
the other effective four-fermion interactions, being only suppressed by the corresponding inverse power of 
the large effective mass scale and eventually by the product of small (Yukawa) couplings. 

2 Although leading order terms in the effective Lagrangian expansion may be actually subdominant in 
definite models. 
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energies in Section |V] Finally, last Section is devoted to conclusions and to comment on 
some necessary ingredients of more fundamental models with RH neutrino interactions. 



II. DIRAC VERSUS MAJORANA NEUTRINOS 

Let us evaluate the suppression factors for the different transitions. The effective Hamil- 
tonian describing the coherent neutrino scattering in a dense medium has the form 

H int = £ £ [ (Za) . . ( ^ rVj) + ^ . . (y.r^)] , (1) 

a i,j 

where T a = I, 75, 7 M , 7^75 and a^ u = 7"], and z a are the complex scalar (25), 
pseudoscalar (zp), vector (zy^), axial- vector (za^) and tensor (zt^u) matrices characterizing 
the background. The sum i,j runs over mass eigenstates. For the moment it is enough to 
know that the Hamiltonian is hermitian, what follows from the relations 

f r a , for a = S, V, A, T, 
r = 70 (r«)t 70 = (2) 
{- T a , for a = P . 

The j i matrix element 

Ti-ij = { v i0^ nt \ v i) 

can be calculated using the plane wave decomposition for the neutrino fields and contracting 
the corresponding creation and annihilation operators with those of the initial and final Dirac 
neutrino, Dirac antineutrino or Majorana neutrino states. We obtain for each case 

= Ej0a)ij UiT a Uj + (Z* a )ji UiT a Uj ] , 
H§ = - E a [(Zah V^ Vl + (Z* a ) tj V^Vt) , (3) 

respectively. Yet it is more convenient to rewrite them only as a function of positive frequency 
spinors using the charge conjugation relation v i = C uf, together with Eq. (T5]) and the 
relations 

T a , for a = S, P, A, 
C{T a ) T C- l = { (4) 
- T a , for a = V, T. 



Then, they read 

Hfj = (Zs)ij UiUj + (Z P )ij tiffsUj + (Zy^ij UiYuj 
+ (Z Afl ) ij ua fl -f 5 u j + (Z Tfll/ )ij Uia^Uj , 

ftfr = (Z^ijUiUj - (Z P )ijUa5Uj - (Z^ijUa^Uj 

(5) 

Hf- = 2(ReZ s )ijUiUj + 2i(ImZp) ij Ui^ 5 u j + 2i(ImZ v ^) ij u i 'y> J -u j 
+2 (ReZ A/Ji )ij ua^uj + 2 i (ImZ Tl J)ij UiO^Uj , 

where the capital Z matrices are constructed from the small z ones. They are hermitian 
Z = z + z 1 for S,V,A,T and antihermitian Z = z — z^ for P terms. As we shall see, Eq. 
(jSJ) is further simplified for relativistic neutrinos. 

A. Relativistic limit 

In the Weyl representation the spinor u\(p u ) for a particle with helicity |, energy E v and 
momentum p u , and normalized in such a way that u*u = 1, has the form 



u\{p v 



2E V 



(6) 



where Xx is the Pauli helicity spinor properly normalized x X — 1- m order to compare 
the spinor products for relativistic neutrinos with different masses, it is convenient to use 
a common energy, which is independent of the particular small neutrino masses involved. 
Indeed, noting that in forward neutrino scattering all momenta are parallel = pi~t and 
following [30], we can write the energy and the momentum of any light neutrino 

where £ is some parameter which depends on the production or detection process, and E 
is the neutrino energy for zero neutrino mass. Using these expressions we can evaluate the 
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(8) 



spinor products entering in Eq. (J5j) 

^(AH(A) = ^ + 0((f) 3 ), 

Ui(X)uj(-X) = 0, 

^(A) 7 S(A) = A^ + 0((f) 3 ), 
« i (A) T 5 « i (-A) = 0, 

^(A) 7 %(A) = r^ + 0((§) 2 ), 

^(-1) 7 %(+1) = -[^(+1) 7 %(-1)]* = m» ^P- + 0((f ) 3 ), 

Mj (A)7^(A) = An^ + 0((§) 2 ), 

^(-1) t m T5U .(+i) = [^( + i) rT5U .(_i)]* = m ^M^ + o((f)3) ) 

^(AKS(A) = *n fc ^p + 0((f) 3 ), 

Ui(X)a kl Uj (X) = A =gp + 0((§) 3 ), 

u<(-l) < 7 0fc u i (+l) = = ?m fc + 0((f ) 2 ), 

^(-1) ( t%(+1) = [« i (+l) < 7 W « i (-l)]* = £ W W + 0((f ) 2 ), 
where = with rf = (sin^cos^, sin^sin^, COS6 1 ) the direction of the neutrino 

momentum, and m M = (0, rrt) with rft = (cos # cos p — i sin y?, cos 9 sin + % cosy?, — sin 6 1 ) 
orthogonal to r? and defined up to an unphysical global phase; whereas e klr is the totally 
antisymmetric tensor with £ 123 = 1. 

Hence, we see that in the relativistic limit for forward neutrino scattering the scalar and 
pseudoscalar terms can be effectively omitted. If in addition we assume that the tensor term 
is negligible (the general case including non-vanishing tensor interactions will be considered 
elsewhere), then not only the corresponding suppressed spin non-flip but the spin flip terms 
can be omitted, and Eq. (jSJ) reduces to 

Wg(A) = (Z v , + XZ A ^ n/> , 

^f(A) = {-Z^ + XZl^.n^ (9) 

Hff(X) = 2 (i ImZ V/l + X ReZ Ay ) i6 n" . 
Thus, the allowed transitions do not flip helicity and depend on the vector and axial-vector 
couplings only. In order to answer the question on the nature of the neutrino mass we must 
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compare the amplitudes TCfj(X = — 1) with TCfj(\ = —1) for particles, and TCR(X = +1) 
with Hfj(X = +1) for antiparticles. 

In the uSM, where neutrinos interact only through LH currents, the vector and axial- 
vector Z matrices are related, Zy^ = —Za^, and the only non-vanishing transitions read 

Hf>(\ = -1) = Wjf (A = -1) = 2(Z v „) ij nf, 

(10) 

HU(\ = +1) = H%(\ = +1) = -2{Z v ^n». 
In any model where neutrinos besides interact through RH currents Zy^ ^ —Za^, and there 
are deviations from the vSM predictions in Eq. (fTOj) which also affect differently the Dirac 
and Majorana neutrino propagation. Indeed, Eq. (Q gives our main input 

ttg(A = -1) - Wjf (A = -1) = (Z* vil + Z\\ n» , 

(11) 

HU(X = +1) - H%(\ = +1) = - (Z V „ + Z Av ) l3 n» . 
Hence, even in the relativistic limit there is in principle the possibility of distinguishing 
between the Dirac and Majorana neutrinos, but only if besides the dominant LH current 
there is some RH piece. Whether this can have practical consequences, it depends on the 
strength of the RH interaction. 

III. LAGRANGIAN AND EFFECTIVE LOW-ENERGY HAMILTONIAN DENSI- 
TIES 

One may adopt a more phenomenological approach at this point and parameterise the 
difference between Dirac and Majorana amplitudes in Eq. (TTTj) . calculating the correspond- 
ing oscillation probabilities afterwards and looking for the largest possible effects. However, 
one must also worry about the physical implications of such a parameterisation and the ex- 
perimental constraints that restrict the different parameters. Then, it is more enlightening 
to start extending the vSM Lagrangian describing the coherent scattering of neutrinos on 
background fermions /, 

Before discussing any specific model, let us introduce an effective Lagrangian with arbitrary 
LH and RH currents and show that the largest effects are associated to new RH neutral 
currents involving the light neutrinos. 
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A. The effective left-right interaction Lagrangian 



The vSM charged current Lagrangian can be generalised to include new LH and RH 
couplings 

e 



C 



cc 



$>v ^(i- 7 5)+^(i+7 5 : 



Kj a w; + h. c . 
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2 v / 2sin6 l H / i>a 

where e£ an d $r are taken to be global factors deviating slightly from their vSM values, 
1 and 0, respectively. This is enough for as we will show, they enter quadratically in 
the effective interaction Hamiltonian. So, |<5#| 2 contributions are negligible because we are 
interested in large (linear) effects. Then, although in contrast with Eq. (fT2|) the RH and 
LH charged currents can have in general different mixing matrices, its explicit form will not 
matter in the RH case. While |£^| 2 terms contribute the same to the propagation of Dirac 
and Majorana neutrinos. Thus, we will take e£ = 1 and U a i to be the 3x3 unitary matrix 
giving the linear combinations of mass eigenstates with well-defined lepton flavour 3l|, [32 ] 



E^ 



One must keep in mind, however, that in definite models the terms neglected can be of the 
same order as those distinguishing between Dirac and Majorana neutrinos below. 
Similarly, the neutral current Lagrangian can be written 

^ " f 75)^ + ^(1 + 75)^, 

75) f}Z» , 



c 



NC 



4 sin 6w cos 6w 



{E 



Vi + 



1 -j 



+ E /V 



7s, 



e N J{l 



(13) 



where as before is a global factor of order 1, its vSM value, to keep track of the order 
of the different terms. The LH piece for light neutrinos can be assumed to be diagonal (and 
e^ u = 1) because it will not help to determine their nature 3 . What is crucial is the general 
form of the RH piece, that we parameterised with a global factor 5r of order 0, the value 
in the uSM, multiplying an otherwise arbitrary hermitian matrix Q^j. The background 
fermions can be taken to be / = e, u, d or e,p, n, with their neutral couplings 



-Nf 



IT) - q f sin 2 6 w + 5 f L , e N R f = -q f sin 2 6 W + 5 f R 



(14) 



In models with neutrino mixing with vSM singlets there are (off-diagonal) terms proportional to this 
mixing and coupling mass eigenstates mainly participating of standard neutrinos and of vSM singlets, 
respectively; as well as quadratic terms correcting the standard couplings. 



s 



where g/ (Tf) is the fermion charge (third component of the weak isospin) and R are 
possible small deviations from their standard values. 

The new RH coupling in Eq. ffT^l) requires the addition of RH neutrinos transforming 
non-trivially under the vSM gauge group, or the Z mixing with an extra gauge boson 
coupling to the new RH neutrinos if these are vSM singlets. Both scenarios have further 
consequences. For example, RH neutrinos in non-trivial SU(2)l representations imply new 
charged leptons which have not been observed, and the Z mixing with an extra gauge boson 
is strongly constrained by processes not involving neutrinos. Definite models must evade 
these constraints. This is not a problem in the Dirac case, but the observed neutrinos 
must mix with the new light neutrinos entering in the RH neutral current if the Majorana 
neutrinos have to feel a different interaction when travelling through matter. Indeed, if 
there are only three light neutrinos feeling the vSM interactions, their nature can not be 
established in oscillation experiments with very energetic neutrinos, although there are new 
sectors beyond the vSM. In order to manifest the Majorana character the light neutrinos 
must participate of new interactions and then of new de gree s of freedom 4 . In either case 



light neutrinos can mixed with new heavy neutrinos |33l . [34J , but this mixing which must 
be rather small is of no relevance to decide about the nature of light neutrinos in oscillation 
experiments. What matters is the mixing with new light neutrinos with RH interactions. 
The experimental limits on such a mixing are not so stringent, because all light neutrinos are 
produced in the standard decay processes and there is no deficit relative to the corresponding 
vSM prediction. On the other hand, in loop processes we have also to sum over all light 
degrees of freedom, and the new contributions are proportional to the (new) light neutrino 
mass differences. 

We could also think in adding new (pseudo) scalar interactions to generate an effective 
low energy Hamiltonian with RH neutrino currents coupled to background fermions, but 
they are also strongly constrained if, as we need, the new interaction couples the observed 
neutrinos to matter fermions. Besides, some (neutrino) mixing with new fermionic degrees 
of freedom and between different spin bosons are still necessary. 



Dirac neutrinos require RH counterparts, but it is not required that these have other interactions. 
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B. The general effective interaction Hamiltonian 



The effective low energy four-fermion Hamiltonian resulting from the former charged and 
neutral interaction Lagrangians has the general form 



n 



eff 



E %E E (^rv,) 



f T a (gja + gjalb) f 



(15) 



f=e,p,n V^2 a =V,A 

where we neglect the scalar and pseudoscalar terms generated by Fierz rearrangement be- 
cause for relativistic neutrinos they vanish or are suppressed by small ^ factors, as discussed 
in Section [Til Thus, we are only left with vector and axial- vector interactions, r a = 7^, 7^75, 
with couplings g l j a and gj a given by 



ij ij 1 ij 

9fV ~ 9fL 1 9fR ? 



-ij -ij 1 -ij 

9fv ~ 9/l ' 9/r j 



(16) 



1 1 

9fA 



-9fL + 9fR 



-11 

9fA 



ii . ii 

-9fL+9fR 



where we introduce LH and RH couplings for later convenience, and 



g% = (A LL + A L % 
g%= {A RR + A R % 



9fL 



_ A LL + A LR^ 



f 

ij ? 



gJ R = (A RR - A RL ) 



(17) 



with 



(A 



LL\f 



-C 12 



where g 



TJ*TJ . X. _L 8 F Nv -Nf r ( ALR\f _ Q Nv p Nf r 

u ei u ej Vfe T 2 L L L V > I 71 Jij ~ 2 L L L R °V 

(A R % = I 5g | 2 U^U ej S fe + § 6%» e N R f Og , (A R % = § 5^e N L f fl 
1. For comparison, within the vSM only 



(18) 



ij 1 



M'z cos 2 e w 



'I I 

9fv 



and 



-9fA 



9fv 



1 1 

-9fA 



(19) 



are different from zero, reproducing the wel 
The form of the Hamiltonian in Eq. (JTB1) 



-known results for LH neutrino interactions. 
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27| is especially useful for our purposes. 



As neutrinos and background fermions / = e,p,n are placed in two separated factors, it 
is straightforward to derive the coherent neutrino scattering in matter summing over all 
fermions and averaging over their properties [32I . This is what is done when going from Eq. 

(dSD to Eq. dU) (aa| 

Wte + i*v%i = - %?<fNf{g%(%) +rn f 9%(i)) , 

( 2 °) 

(z/k + (zYk = (Z/)a = %T. f N f {g%{^)+rn f g%{^)) , 
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where Nf with / = e,p,n stand for the number of fermions / per unit volume, and (-^-) 
and (-gy) are the properly normalised averages over the corresponding fermion distributions 
of the momentum and the spin, respectively. 



C. The effective, coherent neutrino interaction in matter 



The neutrino background coherent scattering is now easy to calculate from the effective 
four- fermion Hamiltonian. Propagation of the light, relativistic neutrinos in matter is gov- 
erned by an evolution equation of the Schrodinger type for each spinor helicity for this is 
conserved (see e.g. [33l. l35||): 

d 



'dt 



(21) 



where \E^( ~]?,t) is the wave function for the neutrino 



antineutrino) eigenstate of mass m*, 



momentum and helicity | = — ~ (+|). As usual [30|, |32|, we assume that all particles have 
the same momentum but different energies Ei = \J ~~ft 2 + rnf . Then, the effective Hamitonian 
in Eq. (T2"TT) reads 



riij 



(P + 



2p' 



(22) 



where TCij describes the coherent neutrino scattering inside matter, thus depending on the 
properties of the medium which is characterised by the fermion contents and the correspond- 
ing polarization, momentum and charge distributions. As it is apparent from Eqs. (Q and 
(I20p . Tiij is also different for Dirac and Majorana neutrinos. In the case of a medium unpo- 
larized, isotropic and neutral (iV e = N p ), the Hamiltonian for Dirac neutrinos with helicity 



is 



Wg (A = -1) = V2G F \N e (g^ L + g%) + N n g^ L 



whereas for Majorana neutrinos of the same helicity it reads 

Wjf (A = -1) = V2 G F \N e (gl 3 L ~ €1 + 9 1 p 3 l ~ 9%) + N n (g^ L - g$) 



(23) 



(24) 



Thus, Dirac and Majorana neutrinos propagate differently in matter in the presence of 
extra RH interactions, g/R^ 0. Their evolution depends on the matter densities Nf and the 
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neutrino couplings to electrons, protons and neutrons gj. These are related to the u and d 
quark couplings in the usual way 

ij r\ ij i ij ij a !j , ij 

9pL — A 9uL ' 9dL > 9pR — A 9uR I 9dR > 

(25) 

ij ij | Q ij ij ij i o ij 

9nL — 9uL i A 9dL i 9nR ~ 9uR ' L 9dR • 

The appropriate Hamitonians for Dirac antineutrinos and Majorana neutrinos of positive 



helicity follow from the well-known relations 33] 



( A ) = ~ [^5(" A )] * > ( A ) = -[< * > ( 26 ) 



which can be also deduced from Eq. ([9]). We should remember that these relations are valid 
in the absence of scalar, pseudoscalar and tensor interactions, for relativistic neutrinos, and 
in the mass eigenstate basis. In the flavour basis these relations are in general not satisfied. 



IV. PROPAGATION OF DIRAC AND MAJORANA NEUTRINOS IN MATTER 

In order to estimate the difference between Dirac and Majorana neutrino oscillations we 
shall consider the simple case discussed in the previous Section of an unpolarized, isotropic 
and neutral medium with constant fermion densities Nf. Then, the neutrino evolution 
equation in Eq. (121]) can be solved analytically. After substracting the non-important 
common diagonal pieces, the Hamiltonians (1231) and ( 1241) can be written in the flavour 
eigenstate basis 

H% (A = -1) = V2G F N e 6 ea 5 eP (27) 

and 

H%(\ = -1) = y/2 G F | N e 6 ea 8 e/3 - N e [ |<5g| 2 6 ea 5 ef3 

+ \ &$T (*F + <*f + + + ft + ft) ng ] (28) 

- N n | 5$" (-1 + 25f + 2S% + 81* + 8%) flg } , 

respectively, with 

n& = E u « fi 3 u h ■ (29) 
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The new RH interactions do not contribute to the Dirac effective Hamiltonian in Eq. (J2] 
which is the same as in the uSM, but they do contribute to the Majorana one in Eq. fl28l) . 
Neglecting quadratic terms in the new small parameters, this reduces to 

= -1) = V2G F N e { S ea S e , + j±\s$>"Q% } • (3°) 

Thus, the difference between Dirac and Majorana neutrino oscillations are triggered by the 
extra term proportional to the neutron to electron density ratio, to 5^ u * and to the Q^a 
matrix elements. Then, to observe any difference, f2ra can not be diagonal. It is also worth 
to emphasize that in contrast with other NP effects, which are quadratic in the (small) 
new parameters, this difference is linear in the (small) strength of the extra RH neutral 
interactions. 

The oscillation probabilities can be found diagonalising the Hamiltonians in Eqs. fl2Tl)and 
( T30l) together with the neutrino (species) dependent kinematical term in Eq. ( |22i) 32 



33| 



W = ±-W {diag (mf)} W/t , (31) 

where W is the diagonalising (unitary) matrix defined by the eigenvectors of 2E v 7i u and 
fhf are the corresponding (real) eigenvalues. We are now ready to calculate the transition 
amplitude from the initial neutrino (production) state \ip(0) > after travelling a distance L 
to some other final neutrino (detection) state \ip(0) >, 

<<p(0)\i(>(L = t)>. 

When considering new neutrino interactions these can also affect their production and 



detection (see e.g. |36|). In general, to calculate the impact of NP on the full process, the 
modification of the initial and final states should be also taken into account. But here we 
are interested in the difference between Dirac and Majorana neutrino propagation, and the 
effects which modify in the same way both type of neutrinos can be ignored. Therefore, we 
can assume that neutrinos are produced and detected in the flavour states \v a > and \vp>, 

respectively. Then the transition amplitude u a — > vp can be written 

~ 2 
777 

A a ^p{L) = < vp{Q)\v a {L = t)>= (W {diag (exp [-i^L])} W^)p a , (32) 



and the transition probability 



29| 



P a ^(L) = 5 aP - 4 ]T R% sin 2 A ab + 8/ 2 J sin A 21 sin A 31 sin A 32 , (33) 



a>b 
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with 



R af3 = Re 



w*w 0a w ab w% 



Fib 



, lip = Im 



w: 2 w P2 w al w t 



pi 



(34) 



and 



A ab = 1.27 



m 2 a - m 2 b )[eV 2 } L[km\ 



(35) 



E v [GeV] 

Using Eq. (1261) we can obtain the transition probability for Dirac antineutrinos and Majorana 
neutrinos with A = +1 from Eq. (13"3"1) with the appropriate replacements 



P^-p{L) = P*^(L; W - W\ G F - -G F ). 



(36) 



V. NUMERICAL RESULTS AND EXPERIMENTAL BOUNDS 



In this Section we quantify how large can the difference between the Dirac and Majorana 
neutrino oscillation probabilities be in the presence of new RH neutral interactions. We 
shall assume that 5r u and the matrices U* and Q R in Eqs. f fT2|) and ([TBI are real, neglecting 
possible CP violating phases because they do not change our results in any substantial way. 
For the calculations shown in the Figures below we use the first order expressions in 5^ u fl R 
for the Dirac and Majorana effective Hamiltonians in Eqs. (I27p and (1301) . respectively. We 



shall take not bigger than 0.01 



"a/3 



37 



(37) 



38l |. and parameterise fl R as follows 

(I 7] if \ 

TJ 1 -x V 
\rj 2 1] 1 — u J 

We have checked that the difference of the Majorana and Dirac neutrino transition proba- 
bilities AP = P AI — P D depends litlle on the diagonal elements (x, for any choice of sign. 
In contrast the dependence on the off-diagonal entry rj is linear 5 . We choose this particular 
f] parameterisation of Q R only for easy comparison of the potential of the different channels 
to discriminate between Dirac and Majorana neutrinos. The oscillation probabilities also 
depend on the medium. We shall concentrate on the neutrino propagation inside Earth, al- 
lowing for variations of the travelling distance L and of the neutrino energy E u . The matter 



5 A quantitative discussion of the effects of flavour diagonal and off-diagonal non-standard Hamiltonian 
contributions to neutrino oscillations can be found in Ref. 26 1. 
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density changes along the Earth radius, and so the electron (proton) and neutron densities in 
Eqs. (1271) and fl30|) . Our calculations are performed for travelling distances L = 1000, 6500 
and 13000 km using along these paths the mean values of the matter densities, which we 
estimate to be equal to 3, 4 and 7 g/cm 3 , respectively. As we are only interested in illustrat- 
ing the main trends, we shall not consider realistic profiles of the Earth matter density. We 
also need to know U. Using the solar, KamLAND and SK+K2K+CHOOZ data, the values 
of the five standard oscillation parameters #12, #13, #23, Sm^ = 5m 2 sol and 5m\ 2 = ^ m ltm can 
be determined. At 95% C.L. [lj 

sin 2 fl 12 = 0.314 (ljftjjj), 

5m 2 21 = 7.92 (1 ±0.09) x 10~ 5 eV 2 , 

sin 2 ^3 = 0.44 (1±££), (38) 

bm\ 2 = 2.4 (l+°; 2 ^) x lO" 3 eV 2 , 

sin 2 9 13 = 0.9 ±g:| x 10~ 2 . 

The numerical results are presented in Figs. 1 — 5. We plot the values of the transition 
probabilities for Dirac and Majorana neutrinos. Their difference AP is of few per cent. Figs. 
1 (a) —1 (f) prove that the transition probabilities also depend on the neutrino flavour. For 

— > v e the effect first increases, and then it seems to weaken again with increasing L. 
But for the — > v T channel the probability difference AP keeps increasing with increasing 
distance L up to the last analyzed value L = 13000 km. As observed in these Figures 
the largest effects manifest in the — > v T transitions when L = 13000 km (approximate 
diameter of the Earth). In this case the largest difference AP between Majorana and Dirac 
neutrino transition probabilities increases with energy, and for rj — 1 reaches its maximum 
at E u pa 57 GeV with a long (experimentally attractive) plateau at higher energies. It varies 
approximately between 0.15 and 0.17, standing for an effect on the value of AP/P D equal 
to 34 % at 50 GeV, increasing up to 84 % at 90 GeV. It can be seen too that for small 
neutrino energies the transition probability difference can be significant as well. However, 
this is of no practical use because of the high frequency of the probability variation, which 
after averaging over the energy bin smooths the signal to zero. Figs. 2 (a) — (b) show the 
77 dependence. The effect weakens with the decrease of rj for it is proportional to S^ u r]. 
Indeed, for 77 = 1/5 the largest difference between the transition probabilities is reached at 
E v pa 53 GeV, with a wide plateau too but with a similar value of AP pa 0.035 (see Fig. 2 
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FIG. 1: Transition probabilities for Majorana and Dirac neutrinos and their difference AP as a 
function of the neutrino energy E v GeV and rj = 1. From left to right and top to bottom: 

(a) U/j, — > u e , L = 1000 km. (b) — ► v T , L = 1000 km. (c) — > v e , L = 6500 km. (d) — ► v T , 
L = 6500 km. (e) -> u e , L = 13000 km. (f) v^^u Tl L = 13000 km. 

(a)). To make it apparent we draw in Fig. 2 (b) the 77 dependance of AP in the vecinity of 
the maximum. Thus, for L = 13000 km and E v = 57 GeV we vary 77 up to 1.0. As expected 
the effect decreases linearly with 77 vanishing for 77 = 0. Finally, in order to visualize the 
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FIG. 2: Left: (a) Transition probabilities for Dirac and Majorana neutrinos and their difference 
AP as a function of the neutrino energy E v in GeV for L = 13000 km and for the process v„ — > v T i 
with 7] = 1/5. Right: (b) Transition probability difference AP for the process — ► v T as a 
function of rj for a neutrino energy E u = 57 GeV and L = 13000 km. 



two-dimensional dependance of AP on L and E u we include Fig. 3. For E u < 20 GeV, AP 
changes very rapidly as a function of the baseline L. For larger energies the variation is 
slower. 

The question arises whether current experimental errors [lO, JT, 12] are not too weak to 



prevent the transition probabilities for Dirac and Majorana neutrinos from overlapping. We 
show in Fig. 4 (a) the — > v T transition probabilities for present values of the oscillation 
parameters (Eq. (1381) ). The bands are obtained varying the oscillation parameters within 
their 95% C.L. limits. Indeed, they are too wide to be able to distinguish between Dirac and 
Majorana neutrinos. However, future experiments will provide more precise measurements. 
Hence, we plot in Fig. 4 (b) the bands for the same central values of the oscillation param- 
eters but assuming that the errors are reduced by a factor of 5. In this case the Majorana 
(upper) and Dirac (lower) bands separate, allowing in principle to distinguish between both 
types of neutrinos. Obviously varying rj (and the form of O) one can make the effect much 
smaller and unobservable. 

This encourages to search for the neutrino character at future experiments. One must be 
aware, however, that not only the experimental errors of the standard oscillation parameters 
must be significantly reduced but they must be determined independently to avoid new 



confusion 



391 ] . A relevant related comment is that the observation of deviations from the SM 
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FIG. 3: Transition probability difference AP as a function of the neutrino energy E v in GeV and 
the distance L in km for the process ZAj — ► z/ T with 77 = 1. The calculations are performed for matter 
densities equal to 3, 4 and 7 [g/cm 3 ] and travelling distances L between — 4500, 5000 — 9000 and 
9500 — 13000 km, respectively. The projections on the L-AP plane are drawn with dots. 



predictions, as those drawn for Majorana neutrinos in the presence of RH neutral currents 
in Fig. HI can have their origin in the NP previously discussed or be a manifestation of the 
inverted character of the neutrino mass hierarchy. The distances and energies relevant in our 
case are similar to those sensitive to atmospheric neutrino oscillations, and then to Sme- 
lts sign, which is unknown at present and may be positive (normal) or negative (inverted 
scheme), can not be fixed in the — > v T channel within the uSM 6 . Therefore in Fig. 
[5] we plot the same transition probabilities as in Fig. [4] (b) but for the inverted scheme. 
The central lines and future bands for Dirac neutrinos, which give the same predictions as 
the uSM, almost coincide for both sign assignments. In contrast, the band for Majorana 



If 6*i3 is sizeable, the sign of 8m} atm will be possibly established in — > v e oscillation experiments 40j . 
even if there are no new interactions. In our case the relevant channel is — > v T because it is where the 
effect of new RH neutral interactions is larger. 
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FIG. 4: Upper and lower bands corresponding to transition probabilities for Dirac and Majorana 
neutrino, respectively, as a function of the neutrino energy E v in GeV for L = 13000 km and the 
process v„ — > v r with 77 = 1. Left: (a) For current 95 % confidence intervals of the oscillation 
parameters [12]. Right: (b) For future 95 % confidence intervals of the oscillation parameters 
(assuming a statistics ~ 25 times larger than nowadays). 

neutrinos moves below (above) the Dirac neutrino band for the inverted (normal) hierarchy, 
reversing the sign of AP but maintaining its size to a large extent. The same displacement 
is obtained changing the sign of rj. This is well-known |4l|. The sign of 5m 2 atrn can be fixed if 
there is other contribution to the interaction potential energy with a well-defined sign, with 
which the standard contribution can interfere. In our case the extra piece is provided by the 
new RH interactions. If in the future any of the two schemes is established, the observation 
of the Majorana character of the light neutrinos in the process in Figs. H] (b) and [5] would 
also allow for the determination of the sign of rj. 

VI. CONCLUSIONS AND MODELS 

Within the vSM oscillations of relativistic neutrinos do not differentiate between Dirac 
and Majorana neutrinos. This is not in general the case beyond the uSM, and Dirac 
and Majorana neutrinos can be differentiated by how they propagate in matter. Indeed, 
extended models with new RH currents or tensor terms do distinguish between both types 
of neutrinos. However, the spin-flip tensor transitions require a polarised medium to manifest 
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FIG. 5: Transition probabilities for Dirac (upper band) and Majorana (lower band) neutrinos in 
the inverse mass scheme (assuming a statistics ~ 25 times larger than nowadays). There is almost 
no difference between the Dirac neutrino bands for the inverted and normal schemes (see Fig. H] 
(b))- 

Moreover, the production and detection of (anti)neutrinos of the wrong helicity are 
strongly suppressed. Hence, we concentrated on the analysis of additional RH currents which 
allow for unsuppressed spin non-flip transitions. This NP does not modify the oscillation 
of Dirac neutrinos in matter for we can assign to them a definite lepton number, but it 
does modify the propagation of Majorana neutrinos for they participate of both types (LH 
and RH) of currents. The largest effect, which is linear in the new small parameters, is 
associated to extra RH neutral currents. However, it is needed a non-trivial flavour structure 
to differentiate Dirac from Majorana neutrino oscillations in matter. Both LH and RH 
neutral current interactions can not be simultaneously diagonal. 

Numerical estimates show that the largest differences between Majorana and Dirac neu- 
trino oscillations manifest in the — > u T channel. This difference between both probabilities 
can be as large as ~ 0.16 for the normal as well as for the inverted scheme, but with opposite 
sign. Unfortunately, this large effect holds for the channel which is more difficult to measure 
and for a very large baseline L ~ 13000 km. Smaller effects are visible for — > u e and other 
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baselines. 

The large effect found is a consequence of the linear, and then unsuppressed, dependence 
of the transition probabilities for Majorana neutrinos on the new RH neutral currents. This, 
which is our main observation, relies on two ingredients for being in practice of some rele- 
vance. The strength of the new flavour violating RH transitions has to be sufficiently large, 
and the precision reached in oscillation experiments sufficiently high. We have not attempted 
a global fit to determine the present experimental limits on the new (flavour changing) cou- 
plings but we have looked for a class of vSM extensions which may accommodate such new 
terms. Let us present a simple example of a class of extended gauge models of possible 



cosmological interest [42j. Consider an extended gauge group with an extra £7(1) jv under 
which all known particles transform trivially, and enlarge the matter content with two new 
fermion singlets under the vSM but with U(1)n charges 2 and —2 (what makes the model 
free of gauge anomalies) and the same lepton number (what forbids direct mass terms in 
the Dirac case). This is completed with two extra Higgs doublets h\ and an extra Higgs 
singlet cf) with U(1)n charges 2, —2 and 1, respectively. The new fermions get their masses 
through their very small Yukawa couplings with the vSM neutrinos, once b! { get non-zero 
vacuum expectation values (vev) v' { , also (in principle) much smaller than the vSM Higgs 
vev < h > = v. Whereas the singlet vev < <fi > = x is (much) larger than v, and gives a 
mass to the new gauge boson. The new RH neutral current term for the light neutrinos is 
proportional to the Z boson mixing with the new gauge boson, which also scales with v\. 
The vSM is practically recovered for < h\ > — 0. Then, the new phenomenology of this 
model depends on those vev. In the absence of Majorana masses lepton number is conserved 
and there are one massless (LH) neutrino and two light Dirac neutrinos. If there are also 
light neutrino masses, we have five light Majorana neutrinos. In the first case the sum in the 
extra RH piece is over the new RH fermions only, while in the second one the sums are over 
all the Majorana neutrinos 7 . One can write down more complicated models with at least 
three massive Dirac neutrinos. Nevertheless, the important question is if the new gauge 
boson and the extra scalars can escape detection and at the same time the LH interactions 
approach the vSM ones, as required by experiment, whereas the new RH piece has flavour 
changing couplings large enough, ~ 100 times smaller than the standard couplings. The 



7 The LH charged current sum also extends to the new light neutrinos. 
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main tension manifests in the Majorana case because the effective RH couplings are the 
product of the Z mixing with the new gauge boson times the square of the mixing of the 
standard neutrinos with the new singlets, and none of them can be negligible. However, the 
usual limits can not be directly applied for the new gauge boson couples weakly (through its 
mixing with the standard Z boson) to charged fermions and the induced flavour violation 
is proportional to the tiny neutrino masses 8 . Obviously, a detailed analysis of the present 
experimental constraints from precise electroweak data |43j is necessary to decide on the 
allowed region of parameters. Such a study is beyond the scope of this paper. 

Our generic conclusion is that we should be aware of the possibility that NP may allow 
to differentiate between Dirac and Majorana neutrinos propagating in matter, and that it is 
worth to search for this difference in the next generation of very precise neutrino oscillation 
experiments. A positive signal would not only indicate the Majorana character of massive 
neutrinos but the existance of new interactions. 
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